Quasiparticle scattering time in superconducting films: from dirty to clean limit 
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We study the quasiparticle energy relaxation processes in superconducting Nb films of different 
thicknesses corresponding to different electron mean free paths in a state far from equilibrium, that is 
the highly dissipative flux-flow state driven up to the instability point. From the measured current- 
voltage curves we derive the vortex critical velocity v* for several temperatures. From the v*(T) 
values, the quasiparticle energy relaxation time r £ is evaluated within the Larkin-Ovchinnikov model 
and numerical calculations of the quasiparticle energy relaxation rates are carried out to support 
the experimental findings. Besides the expected constant behavior of r e (T) for the dirty samples, 
we observe a strong temperature dependence of the quasiparticle energy relaxation time in the 
clean samples. This feature is associated with the increasing contribution from the electron-phonon 
scattering process. 
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I. INTRODUCTION 



The quasiparticle energy relaxation time is one of 
the fundamental properties of superconductors which 
are essentially non-equilibrium in nature. Despite the 
great interest towards the comprehension of such mecha- 
nism, there are still few experimental techniques avail- 
able to measure quasiparticle lifetimes and most of 
these are performed when the superconductor is slightly 
far from equilibrium as for example the quasiparticle 
injection^ the photon irradiation^ or the measure of 
the quasiparticle-lifetime broadening in superconducting 
tunnel junctions^ 

The study of the electronic flux-flow instability provides 
a way to extract information about the quasiparticle en- 
ergy relaxation from an intrinsic bulk phenomena of the 
material when the quasiparticle distribution is far from 
equilibrium4rii Indeed, the electronic flux-flow instabil- 
ity is a phenomenon associated with high velocity vortex 
motion, a dynamic regime characterized by a continuous 
diffusion of quasiparticles from the vortex core to the sur- 
rounding superconducting medium and by a re-trapping 
of quasiparticles into the vortex coreZ 
Talking about the mechanisms which lead to the elec- 
tronic flux-flow instability, we have to distinguish two 
limits: the high-temperature limit (T w T c ) and the 
low-temperature limit (T <C T c ). Close to the critical 
temperature the Larkin-Ovchinnikov (LO) theory can be 
applied^ where the instability is associated to a shrinking 
of the moving vortex and to a drastic drop in the viscos- 
ity of superconducting medium due to the diffusion of 
quasiparticles from the vortex core to the surrounding. 9 
On the contrary, in the low-temperature limit there are 
three effects, i.e. the raising of the electronic temper- 
ature, the creation of additional quasiparticles and the 



diminishing of the superconducting gap A, which lead 
to an expansion of the vortex*^ This expansion has the 
effect to reduce the viscous drag because of the softening 
of gradients of the vortex profile ^ 

Although the two temperature regimes are different, the 
flux-flow instability signature in the current- voltage char- 
acteristics is the same: an S-shape if voltage driven mea- 
surements are performed, or a sudden jump in voltage 
in the current driven ones. Moreover the instability pa- 
rameters are directly related to the quasiparticle energy 
relaxation time r e . In particular, in the LO framework it 
is: 



(1) 



where D is the quasiparticle diffusion coefficient, ((x) is 
the Riemann function and v* is the vortex critical veloc- 
ity. 

Finally, it is known that the quasiparticle energy relax- 
ation process is ruled by different mechanisms in the dirty 
and clean limits In the dirty limit, the energy relaxation 
occurs within the vortex core, while in the opposite sit- 
uation (clean limit) quasiparticles can relax the energy 
absorbed by the electric field only after experiencing a se- 
ries of Andreev reflections and covering a total distance 
within the core equivalent to several coherence lengths 
Vs. 

In this work we analyze the temperature dependence of 
the quasiparticle energy relaxation time as deduced by 
the flux-flow instability in a fixed magnetic field and in 
the range of temperature from « 0.5T C to ~ T c on a 
set of Nb superconducting strips of increasing thickness 
to grasp the variation from the dirty to the clean limit. 
To support the experimental results a combined analy- 
sis based on the LO theory and the well known model 
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of quasiparticle energy relaxation developed by Kaplan 
et alJ^ is carried out in the whole range of tempera- 
ture to describe the crossover from dirty to clean limit. 
We observe that the temperature range of our experi- 
ment (4.2 10 K) is well above the temperature value 
T x = r 2 ksTp/eF, where To is the Debye temperature, 
cf is the Fermi energy and r is the phonon reflection coef- 
ficient at the film-substrate interface arising from acous- 
tic mismatch and it is always r < l. 6 This temperature 
is the threshold above which the electron-electron scat- 
tering time is higher than the electron-phonon one and, 
consequently, the LO approach can be retained. We esti- 
mate from known parameters^ that for Nb this crossover 
temperature is T x = l.lr 2 K. 

We find that while the total quasiparticle relaxation time 
decreases toward the dirty limit, the electron-phonon 
scattering contribution increases with respect to the re- 
combination one. 



II. QUASIPARTICLE ENERGY RELAXATION 
TIME 

Usually, quasiparticle energy relaxation rate t~ is 
considered as the result of the combination of differ- 
ent contributions from three different classes of pro- 



cesses: the electron-phonon scattering, the recombina- 
tion to Cooper pairs and the inelastic part of electron- 
impurity scattering. The first kind of processes involves 
the changing of the quasiparticle excitation energy by 
means of interaction with phonon, while the second leads 
to the emission of a phonon by the two recombining par- 
ticles to form bound Cooper pairsJ^ In the limit of low 
concentration of non-magnetic impurities the electron- 
impurity scattering is usually neglected, indeed the prop- 
erties of an s-wave superconductor are quite insensitive 
to this scattering process which plays a role only at very 
low temperature^ 

The quasiparticlcs lifetimes due to the first and sec- 
ond type of processes described above are related to 
phonon density of states, F(SY) (f2 being the phonon fre- 
quency), and the matrix element of the electron-phonon 
interaction, a 2 (SI), within the strong coupling Eliashberg 
formulatio n 15 ' 16 as discussed in Ref. [l2]- Within this 
approach, the inverse of the lifetime, t~ (cj) (ui being 
the quasiparticle frequency), is equal to twice the decay 
rate T(w), which can be calculated by the single-particle 
Green's function^ 2 , or by means of the golden rule, i.e., by 
calculating the probability that the electron in the state 
with wavevector p and frequency oj will emit a phonon 
with wavevector q and frequency f2, and its explicit ex- 
pression is: 



^' M-smJr dsW < n > F(ii)Re (tf^M - A) m ~ "> + " (n)l + 



where f(io) and n{uj) are the Fermi and Bose distri- 
bution functions, respectively, while the spectral weight 
a 2 (Q)F(Q) is given by averaging the square of the dressed 
electron-phonon matrix element over the Fermi surface, 
Zi(0) is the renormalization parameter and A is the su- 
perconducting gap. 

The first and third term in Eq. [^represent the scattering 
processes with the emission and absorption of a phonon, 
respectively; they are denoted with r s . The second term 
describes the recombination of two quasiparticles to form 
a pair and whose excess energy is emitted as a phonon; 
it is denoted with r r . For ui ~ kT c and T <T C only the 
low-frequency part of a 2 F contributes where it behaves 
quadratically as w Ml 2 , with 10 3 & equal to 4 meV~ 2 for 
the Nb.— The above expressions of r s and t t will be 
employed in the comparison with the quasiparticle relax- 
ation times deduced from the vortex critical velocities. 
Despite the fact that both electron-phonon scattering 



and recombination processes transfer directly the excess 
quasiparticle energy to the crystal lattice and (eventu- 
ally) to the heat bath, the temperature dependence is 
distinctly different for both processes. Indeed, the scat- 
tering lifetime r s associated to the electron-phonon scat- 
tering increases as the temperature is lowered, with a 
power law of the type r s cx T~ n with n w 3, owing to 
the decrease of the phonon population. At the gap edge 
uj = A(T) the quasiparticle cannot emit a phonon and 
scatter because it is in the lowest energy state. For quasi- 
particles with energies u > A(T), spontaneous emission 
of a phonon sets a limit to the scattering lifetime r a , thus 
t s becomes almost temperature independent^ 2 - 
The recombination lifetime, instead, must reflect the ex- 
ponential temperature dependence of the quasiparticle 
population e A / kT . Another peculiarity is that for a 
quasiparticle at the gap edge t t goes through a mini- 
mum value for temperature close to T c , while for w > A 
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TABLE I. Physical parameters of the Nb films. 





Sample 


w [/jm] 


d [nm] 


Tc [K] 


I [nm] 


13 


dirtiest 


NbA02 


40 


30 


7.6 


1.6 


1.7 


I 


NbA3 


50 


60 


8.1 


2.3 


2.1 


I 


NbC6 


50 


135 


8.7 


2.3 


2.3 


cleanest 


NbD8 


50 


150 


9.2 


3.7 


2.6 



there is no minimum ^ 

Let us note that when T approaches T c , the limit- 
ing values of t s and r r for ui — A(T C ) = are 
equal due to particle-hole symmetry and their value is 
T s{r) (0,T c )/T = 4.20. Here r = Z 1 (0)h/2nb(kT c ) 3 is a 
characteristic time which is used as a time unit to express 
t s and T r in an universal form^ 

Finally, we must consider that there are situations in 
which the branches of the quasiparticle excitation curve 
corresponding to quasiparticle wavevectors less than and 
greater than the Fermi wavevector are not equally occu- 
pied. In these situations, one can distinguish among the 
processes we are considering those which contribute to re- 
laxation of the branch population imbalance. The associ- 
ated lifetime is called the branch-mixing time tq . There 
is a relation between tq, t t and r s : Tq 1 < | (tjT 1 + t" 1 ) 
(where the equality holds at the gap edge ui = A).— 



III. EXPERIMENTAL DETAILS AND RESULTS 

The samples investigated are Nb strips patterned by 
a standard UV photo-lithography technique on Nb films 
deposited on Si(100) substrates in a UHV dc diode mag- 
netron sputtering system. The film thickness d was con- 
trolled via a quartz crystal monitor calibrated by low- 
angle reflectivity measurements. The Nb strips width w 
is 40 fim for the sample named NbA02 and 50 /jm for 
samples NbA3, NbC6 and NbDStf. The distance L be- 
tween the voltage tips for all the samples is 2 mm long. 
A typical value of the critical current density, obtained 
by the standard 10 /iV/cm criterion, at 4.2 K and zero 
magnetic field is J c = 2 • 10 6 A/cm 2 . 
The characteristic parameters of the samples are sum- 
marized in Table HI The electron mean free path I is 
estimated via the relation I — 1 ' — 



v F 



vfJPn V e J 

2.73 • 10 5 m/s is the Fermi velocity^ 



— where 
7 = 7- 



10 2 J/K 2 m 3 is the Nb electronic specific heat coefficient 20 
and the resistivity pm is estimated at T = 10 K. The 
mean value of the residual resistance ratio is around 
(3 = R(300K)/R(10K) w 2, while it increases as the 
thickness increases (see Table HI thus the samples are get- 
ting progressively clean with increasing thickness. 
As it is well known, when £o <C I we are dealing with the 
clean limit, while for £ 3> I the sample is in the dirty 
limit. Thus, comparing the values of I in Table|I]with the 
value £o ss 7 nm of the coherence length as obtained by 
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o.s 
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FIG. 1. Experimental v*(t) curves at B = 0.01 T of sam- 
ples NbA02 (orange triangles), NbA3 (red circles) and NbD8 
(blue squares). The solid lines are the resulting fitting curves 
obtained by Eq. assuming for r e the form in Eq. [H where 
t s and r e i are the fitting parameters. 



H C 2 (T) measurements and assumed to be independent of 
djSl we can argue that the whole set of samples covers 
the crossover region between the two limits. The reduc- 
tion of the mean free path with decreasing thickness of 
the samples may be due to the poor crystalline quality of 
Nb film in proximity of the interface with the substrate. 

In order to analyze the mean free path influence on 
the quasiparticle energy relaxation time, we use as a tool 
the flux-flow instability phenomenon which is observed 
in the V{I) curves of the samples for a wide tempera- 
ture range j 17 ' 22 In particular, in Fig. [T] are shown the 
v*(t) (with t = T/T c ) for the three samples of increasing 
I: NbA02, NbA3 and NbD8. The curves are extracted 
from current driven V(I) measurements performed in 
a fixed magnetic field B = 0.01 T with a pulsed cur- 
rent 4-probe technique in order to minimize self-heating 
effects.— The v* values are evaluated via the formula 
v* =V*/(BL), where V* is the voltage value at the in- 
stability point. For the dirty samples NbA3 and NbA02 
the v* (t) shows a monotonous decrease of the vortex crit- 
ical velocity as the temperature is raised, in agreement 
with the LO expectation^ On the contrary, the v* (t) for 
the cleanest sample NbD8 presents the evident existence 
of a maximum before the predicted decrease. This un- 
usual behavior can be ascribed to the clean nature of the 
considered sample, as it is explained in Section IV. 
Let us note that, even if the temperature range is not 
so close to T c , the use of the LO approach is further 
sustained by the evaluation of the value of the Bezug- 
lyj and Shklovskij parameter Indeed, the Bezug- 

lyj and Shklovskij (BS) approach extends the LO model 
taking into account quasi-particle heating due to the fi- 
nite heat removal rate of power dissipated in the sample. 
In the BS approach a macroscopic parameter Bt is de- 
rived from a microscopic analysis of the heat removal as 
a result of the mutual phonon exchange rather than in- 
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FIG. 2. Dissipated power density p* as a function of the 
applied magnetic field B at T = 4.2K for the same three 
samples of Fig. [1] The solid lines are guide for the eyes. 



terface properties between the film and the substrate 24 . 
The Bt parameter separates the region where non ther- 
mal (B -C Bt) or pure heating mechanisms (B 3> Bt) 
of the instability dominates. In our measurements Bt 
at 4.2 K is 0.24 T,— i.e. a value well above the applied 
magnetic field B = 0.01 T. Thus, the quasiparticle over- 
heating, which is effective only when the magnetic field 
value exceed Bt, is negligible. 

Moreover, the magnetic field dependence of the dissi- 
pated power density p* — J* ■ E* is considered (see Fig. 
[3]), where J* is the current density and E* is the elec- 
tric field value at which the instability occurs. There is a 
monotonous increase of p* with the increasing field which 
indicates that the self- heating effects are negligible^ 



IV. DISCUSSION 



spatial homogeneity of the non-equilibrium quasiparticle 
distribution between the vortices in the range of low mag- 
netic fields^* 2 ! 

From Fig. [3] it is clear that for all three dirty samples 
(namely, NbA02, NbA3 and NbC6) an almost tempera- 
ture independent r n (t) is observed. The slight increase 
related to samples NbA02 and NbA3 for t > 0.9 will be 
discussed in the part B of this Section. 
On the contrary, in Fig. [T]a strong non-monotonic behav- 
ior of the v* (t) related to cleanest sample NbD8 can be 
clearly seen. This behavior is associated with a temper- 
ature dependent r e (see Fig. [3]). Indeed, in this limit the 
contribution of the recombination processes is not neg- 
ligible and it is known that these processes are strongly 
influenced by the temperature (see Section ITT| . 
In order to estimate the contribution of the recombina- 
tion processes we performed a non-linear fitting proce- 
dure on the data shown in Fig. [T]on the basis of Eq. |3jby 
assuming the quasiparticle energy relaxation rate r" 1 as 
the sum of two contributions: one is a temperature inde- 
pendent term t~ , which takes into account the scatter- 
ing processes, and the other one is an exponentially tem- 



perature dependent term r r 



i'Exp 



-1.76 to ■ 



which takes into account the recombination processes: 



1 

TV 



1 



1 

Tel 



-Exp 



-1.76 to 



t 



(4) 



where to is a numerical parameter. The assumption to 
consider the scattering contribution as independent from 
temperature follows from the observation that when the 
instability is triggered, the quasiparticle energy is much 
higher than the energy gap, as we demonstrate below. 
In Table |TT] are summarized the results of the fitting pro- 
cedure. The values of the two fitting parameters r e z and 
t s were obtained by fixing the parameter to equal to 1.5, 
which can be interpreted in the framework of recombina- 



A. On experimental data 



The experimental data presented in Fig. [T] related to 
the dirty samples NbA02 and NbA3 fairly agree with the 
LO prediction on a monotonous decrease of the vortex 
critical velocity as the temperature approaches T c . This 
behavior is associated with a temperature independent 
quasiparticle energy relaxation rate t" 1 4 This prediction 
is confirmed by the evaluated r n (t) = T e (t)/r e (0.5) curves 
shown in Fig. [3J where the T e (t) are estimated from the 
experimental v*(t) curves inverting the relation: 



(1 -t) 1/2 D [14C(3) 



1/2 



7T T e 



«0 



(3) 



with arj 



D 



\vfI- 



1/2 



the inter-vortex distance and 



( JL*o 
VV3 b 

Eq. UJ is a derived expression from the 



LO which takes into account the necessary condition of 




FIG. 3. Evaluated r„(t) curves by Eq. CD at B = 0.01 T for 
samples NbA02 (orange triangles), NbA3 (red circles), NbC6 
(green stars) and NbD8 (blue squares). Solid line are guide 
for eyes. 
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TABLE II. Estimated quasiparticle energy relaxation times. 





Sample 


Tel [ns] 


t 3 [ns] 


T e (t = 0.5) [ns] 


t [ns] 


p*(t = 0.5) [10 12 W/m 3 ] 


dirtiest 


NbA02 


« 5 


0.33 ±0.01 


0.32 


0.26 


0.58 


I 


NbA3 


1.60 ±0.90 


0.42 ±0.01 


0.42 


0.22 


0.54 


cleanest 


NbD8 


0.22 ±0.01 


0.71 ±0.01 


0.66 


0.15 


0.35 



tion processing involving more than two quasiparticles.-i 
If we compare for each sample the t s and r e i values, we 
recognize that for the dirty samples (namely NbA02 and 
NbA3) the faster process is the scattering one, thus it 
is the one which dominates the quasiparticle energy re- 
laxation and it results r e w r s . On the contrary, for 
the cleanest sample (namely NbD8) it is the recombina- 
tion process which dominates over the scattering one. As 
shown in Fig. [TJ the exponential dependence on tempera- 
ture of the recombination processes determines the pecu- 
liar non monotonic behavior of the vortex critical velocity 
in this sample. These results point out that the contribu- 
tion from the recombination process becomes more sig- 
nificant as the sample becomes more clean. 
Theoretical and experimental studies on the temperature 
dependence of the lifetime broadening factor T of the 
tunneling conductance curves gave also evidence of an 
electron-phonon scattering rate larger than the recombi- 
nation one in disordered superconductors^ - — Further- 
more a power law dependence of the scattering rate on 
temperature (~ T 3 ) has been deduced from these tunnel- 
ing experiments, whereas an almost constant behavior of 
t s can be inferred from our measurements. Neverthe- 
less, in agreement with a previous study on the flux flow 
properties of clean and dirty samples,— here we find the 
experimental evidence of a crossover between two dif- 
ferent behaviors: for dirty samples the electron-phonon 
scattering dominates, while for the clean samples the re- 
combination to Cooper pairs plays the significant role, 
and the change in the balance between the two contribu- 
tions is continuous. 

Moving from the dirty to the clean limit r e ; decreases, 
while r s increases as well as the total quasiparticle en- 
ergy relaxation time t £ (evaluated via Eq. as it can 
be seen in Table [TTJ The decrease of the scattering time 
t s in the thinner samples and its temperature behavior 
may be ascribed to an additional relaxation channel due 
to the scattering by impurities and/or to the high en- 
ergy quasiparticles involved in the conventional electron- 
phonon interaction. Indeed, whereas in first approxima- 
tion the energy gap and the thermodynamic properties 
of superconductors remain unchanged in the presence of 
a low concentration of impurities,— our results may indi- 
cate that the situation is more complex and the inelastic 
part of the impurity scattering may play a role. 
On the other hand, in the limit of high quasiparticle ex- 
citation energy a drastic decrease and a nearly constant 
temperature dependence of the electron-phonon scatter- 
ing lifetime is predicted^ In Fig. 0] there is a plot of the 
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ci : .55 \ 




ca= 1.69A 




co=1.78A 





t 

FIG. 4. Scattering time r s in units of to as a function of 
the reduced temperature t = T/T c for different values of the 
quasiparticle energy: oj = 1.55A (dotted line), 1.69A (dashed 
line) and 1.78 A (solid line). 



scattering time T s (t)/ro for different values of the ratio 
t s (t — > 0) /to (with t s (t — > 0) assumed to be almost equal 
to the estimated t s in Table [TTJ) , which correspond to 
different quasiparticle energies above the gap. As shown, 
the more the sample is dirty, the more the scattering time 
is lowered. 

To support this scenario, an indication of the energy scale 
of quasiparticle excitations involved at the vortex insta- 
bility point comes from the density of electric power p* 
for the three investigated samples (see Table [TTJ) . The p* 
values increase going towards the dirty sample and this is 
indicative of the presence of an high energy quasiparticle 
population. 

This is further confirmed by comparing the experimental 
data of v*(t) with the curve obtained by Eq. fusing the 
expression of the scattering time t s resulting from Eq. [21 
as shown in Fig. [5] To compare the two curves, one has 
to determine the quasiparticle energy (above the gap) by 
the knowledge of the ratio T s (t — > 0)/to. In particular 
the solid line in Fig. [5] includes a scattering lifetime t s 
for quasiparticles with excitation energy uj = 1.78A, well 
above the gap edge. 

Thus, it appears that the electron-phonon scattering con- 
tribution in dirty samples prevails on the recombination 
one, and its weaker temperature behavior deduced from 
our data seems to be related to the quasiparticles from 
different energy scales involved in the electronic flux-flow 
instability mechanism. 
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FIG. 5. Comparison between the experimental v*(t) curve of 
sample NbA02 (dots) and the one obtained by Eq. with 
a value of t s corresponding to quasiparticles with excitation 
energy to = 1.78A (solid line). 




I 

FIG. 6. Branch-mixing time tq in units of to as a function of 
the reduced temperature t = T/T c for oj — A(0) (solid line) 
and oj = A(T) (dashed line). 



B. Other contributions to the quasiparticles 
scattering rates: the branch mixing 

As discussed previously, the observed T e (t) behavior 
is in good agreement for the dirty and clean samples 
with literature However, there is a slight discrepancy 
in the T e (t) curves of the dirty samples NbA02 and 
NbA3, which show a clear up bended curvature as the 
temperature approaches T c (see Fig. [3]). This can be 
interpreted as a signature of the presence of additional 
scattering processes associated to a highly unbalanced 
distribution of quasiparticles in the levels below and 
above the Fermi level, the so called branch-mixing 
process. 

In fact, while in thermal equilibrium the branches of 
quasiparticles with wavevector below and above the 
Fermi surface are equally populated, in the presence of 
a current flow across an interface or due to the injection 
of quasiparticles through tunnel barriers, a branch 
population imbalance can occur. The contribution to 
the quasiparticle decay rate due to the branch mixing 
can be calculated by the golden rule or by keeping 
extra terms in the Nambu self-energy^ 2 In particular, 
when the quadratic approximation is used for the factor 
a(Q)F(n) and the BCS relation 2A = 3.52fcT c is con- 
sidered, a useful approximate form of tq for oj 3> A(T) is: 

^T)^fe) G^) coth (2fcr)- (5) 

As shown tq diverges as A _1 (T) for T — > T c . For other 
ranges of the excitation energy oj the following relation 
holds, Tq 1 < ^(t^ 1 + TjT 1 ), where the equality holds at 
the gap edge oj = A(w). This leads to a limiting value 
of Tq 1 (A(T)) at T = T c given by 0.73. A plot of t q is 
reported in Fig. |B] for oj = A(T) and oj = A(0). One 



notices that at the band edge a minimum develops close 
to T c . 

We can suppose that the mechanism responsible for the 
unbalanced distribution of quasiparticles in the levels be- 
low and above the Fermi level leading to the branch mix- 
ing is the same mechanism responsible for the instability 
in the dirty samples. 



V. CONCLUSIONS 

Summarizing, we have investigated the quasiparticle 
energy relaxation processes from dirty to clean supercon- 
ducting films by studying the temperature dependence of 
the electronic flux-flow instabilities. We have found that 
the vortex critical velocity behavior as a function of the 
temperature v*(t) and the resulting quasiparticle energy 
relaxation time r e (t) are strongly modified when the sam- 
ple changes from the dirty to the clean limit. From our 
findings we can argue that the vortex instability mecha- 
nism involves different quasiparticle energy scales in Nb 
films with different electron mean free paths, giving rise 
to a more important contribution of the electron-phonon 
scattering on the energy relaxation time in dirty samples. 
It appears that the fraction of the quasiparticles experi- 
encing scattering is eliminated from the pairing process. 
These conclusions are supported by numerical calcula- 
tion based on the model developed by Kaplan et al. for 
the quasiparticle energy relaxation processes^ Indeed, 
the simulations show how the v*(t) related to the dirty 
samples can be reconstructed considering only the con- 
tribution from the quasiparticle scattering, while for the 
clean limit it is necessary to consider also the contribu- 
tion by the quasiparticle recombination. 
Moreover, our investigation enlightens a deviation of T e (t) 
from the usually constant behavior of the extreme dirty 
samples. Indeed, the r e (t) curves show a bend up curva- 
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ture for T sa T c . This unexpected behavior is associated 
to a non-negligible contribution to the quasiparticle en- 
ergy relaxation from the branch-mixing processes. It is 
supposed that the cause of the unbalancing distribution 
of quasiparticles in the levels close to the Fermi level is 
the instability mechanism. 

Recently, we became aware of overlapping, simultaneous 
workpi in which the energy relaxation time in low pinning 
molybdenum-germanium films has been measured by the 
vortex instability tool. We hope that further experimen- 
tal and theoretical work will be carried out in order to 
deeper clarify this subject. 
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